Localization of acoustic waves in ID random liquid media 
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We study acoustic propagation in one dimensional water ducts containing many air-filled blocks. 
The acoustic band structures for the periodic arrangements of the blocks is calculated, whereas the 
transmission for various random configurations of the blocks is computed by the transfer matrix 
method. The results show that while all waves are localized for any given amount of disorders, there 
is no genuine scaling behavior for the system. The results also reveal a distinct collective behaviour 
for localized waves, a feature useful for distinguishing the localization from the residual absorption 
effect. 

PACS numbers: 43.20., 71.55J, 03.40K 



The fact that the electronic localization in disordered 
systems Q is of wave nature has led to suggestion that 
classical waves could be similarly localized in random sys- 
tems. The effort in searching for localization of classi- 
cal waves such as acoustic and electro-magnetic waves is 
tremendous. It has drawn intensive attentions from both 
theorists (| and experimentalists |i"l"|Jp| |. 

Localization of waves in one dimensional (ID) systems 
has attracted particular interest from scientists because 
in higher dimensions the interaction between waves and 
scatterers is so complicated that the theoretical compu- 
tation is rather involved and most solutions require a 
series of approximations which are not always justified, 
making it difficult to relate theoretical predictions to ex- 
perimental observations. Yet wave localization in one 
dimension (ID) poses a more manageable problem which 
can be tackled in an exact manner by the transfer ma- 
trix method. Moreover, results from ID can provide in- 
sight to the problem of wave localization in general and 
are suitable for testing various ideas. Indeed, over the 
past decades considerable progress has been made in un- 
derstanding the localization behavior in ID disordered 
systems [[13|. However, a number of important issues re- 
mained untouched. These issues include, for example, 
how waves are localized inside the media and whether 
there is a distinct feature for wave localization which 
would allow to differentiate the localization from residual 
absorption effect without ambiguity @ jl^] ■ Results from 
the statistical analysis of the scaling behavior in ID ran- 
dom media is not conclusive. A further question could be 
whether the localized state is a phase state which would 
accommodate a more systematic interpretation from the 
view of a symmetry breaking and collective behavior, in 
analogy to phase states such as superconductivity. This 
Letter attempts to provide insight to these questions. 

Here we study the problem of acoustic wave propaga- 
tion in one dimensional water ducts containing many air 
blocks either regularly or randomly but on average reg- 
ularly distributed inside the ducts. The frequency band 
structures and wave transmission are computed numeri- 



cally. We show that while our results affirm the previous 
claim that all waves are localized inside an ID medium 
with any amount of disorder, there are, however, a few 
distinctive features in our results. Among them, in con- 
trast to optical cases EtI] , there is no universal scaling be- 
havior in the present system. In addition, when waves are 
localized, a collective behaviour of the system emerges. 

Assume that N air blocks of identical thickness a are 
placed regularly or randomly in a water duct with length 
L measured from the left boundary of the duct (LB). The 
air fraction is f3 = Na/L, the average distance between 
two adjacent water layers is (d) = L/N = a/ (3, and the 
average thickness of water layers is (b) — EjLi bj)/N. 
The degree of randomness for the system is controlled 
by a parameter A in such a way that the thickness of 
the j-th water layer is bj = (b)(l + Sj) with Sj being a 
random number within the interval [—A, A]; the regular 
case corresponds to A = 0. An acoustic source placed at 
LB generates monochromatic waves with an oscillation 
v(t) — ve~ lu}t . Transmitted waves propagate through 
the N air layers and travel to the right infinity. In order 
to avoid unnecessary confusion, possible effects from sur- 
face tension, viscosity or any absorption are neglected. 
For convenience, we use the dimensionless quantity k(b) 
to measure the frequency, where k — oj/c is the wave 
number and c is the sound speed in water. Similarly, k g 
and c g represent the wave number and sound speed in 
the air blocks respectively. 

The wave propagation in such a system can be solved 
using the transfer matrix method [Q . After dropping out 
the time factor e~ luJt , the wave propagation is governed 
by two equations jL6). The first is the Helmholtz equa- 
tion, 

p'n(x) + k 2 rnPm {x) = 0, (1) 

in which p m (x) is the pressure field, and the subscript 
m refers to the medium that can be either water or 
air, depending on where x is located. Within any layer, 
Eq. (Q) warrants two solutions: A m e lkmX represents the 
wave transmitted away from the source to the right and 
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B m e~ lkmX the wave reflected towards the source. The 
total wave is therefore p m (x) = A m e lkmX + B m e~ lkmX . 
The second equation relates the oscillation velocity and 
the pressure field, 
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where p m refers to the equilibrium mass density of 
medium m. 

The coefficients A m and B m in any two adjacent layers 
are connected by a transfer matrix. By invoking the con- 
dition that the pressure and velocity fields are continu- 
ous across the interfaces separating water and air and the 
condition that there is no reflected wave to the right end 
of the system, the matrix elements for the transfer ma- 
trices linking all interfaces can be deduced, and waves in 
any particular block is therefore completely determined. 
The ratio between the amplitude of the outgoing wave at 
the right boundary and that of the transmitted wave at 
the source defines the transmission rate. 

Before going further, a general discussion on wave 
propagation is in place. When waves propagate through 
media alternated with different material compositions, 
multiple scattering of waves is established by an infi- 
nite recursive pattern of rescattering. In terms of wave 
fields, the energy flow in the system is calculated from 
J - Re[i(p*(x)d x p(x)]. Writing p(x) = \A(x)\e i6 W with 
\A\ and 9 being the amplitude and phase respectively, the 
energy flow becomes J ~ \A\ 2 d x 9. Obviously, the energy 
flow will come to a complete halt and the waves could be 
localized in space when phase 9 is constant and \A\ does 
not equal zero. 

First consider the case with the periodic placement of 
air blocks (with spatial period d, water layer thickness b) . 
According to Bloch theorem, and the wave field p can be 
written as 



p[x) — A(x)e 



iKx 
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where A(x) is a periodic function satisfying A(x + d) = 
A(x), and K is the Bloch wave number represented in 
the dispersion relation 



cos K d = cos k g a cos kb 



cosh2?7 sinfcga sinfcfr. (4) 



Here k g = k/h and r\ = lnq with q 2 = gh and g = p g /p, 

h = Cg/C. 

The band structures for four air-fractions are shown in 
Fig. 0. The ranges covered by the dispersion curves refer 
to the passing bands, while the areas sandwiched by any 
two curves to the complete band gaps. Within the gap, 
waves are evanescent and cannot propagate. Since the 
air blocks are strong scatterers due to the large contrast 
in acoustic impedance between air and water, we see that 
even a small air-fraction can lead to wide band gaps as 



seen m Fig. 0(a). With increasing /?, the passing bands 
are narrowed and become streaks. When the air-fraction 
is reduced, however, the band gaps gradually disappear, 
as shown by Fig. 0(d). 



(a) b=io- 



(b) p=10" 



-0.5 0.5 1 

(c) P=10" 5 





FIG. 1. Band structures of an infinite periodic system for 
four different air fractions (/3). In the plots, K is the Bloch 
wave number and d is the lattice spacing. 

Wave propagation properties can be significantly af- 
fected by varying air-fraction or adding randomness. 
Fig. 0(a) presents the typical results of the transmission 
rate as a function of k(b) for various (3 at a given random- 
ness. At frequencies for which the wavelength is smaller 
than the averaged distance between air blocks, the trans- 
mission is significantly reduced by increasing air fraction. 



(a) Transmission vs Frequency (A=0.3 ) 




(b) Transmission vs Frequency (p=1e-4 ) 




FIG. 2. Transmission versus k(b) /it for various air fractions 
at A = 0.3 (a) and different disorders (b). The number of the 
air-blocks is 100. 

Fig. 0(b) illustrates the effect of the randomness A 
on transmission for a given air-fraction. For compari- 
son, the transmission in the corresponding regular array 
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(A = 0) is also plotted. The gaps are located between 
k(b)/it = 0.46 and 1, 1.21 and 2, 2.12 and 3, and so on. 
We find that for frequencies located inside the band gaps 
of the corresponding regular array, the disorder-induced 
localization effect competes yet reduces the band gap ef- 
fect. To characterize wave localization in this case, both 
the band gap and the disorder effects should be consid- 
ered, supporting the two parameter scaling theory Jl7| . 
However, increasing disorder tends to smear out the band 
structures. When exceeding a certain amount, the effect 
from the disorder suppresses the band gap effect com- 
pletely, and there is no distinction between the localiza- 
tion at frequencies within and outside the band gaps. 

Fig. || shows that the localization behaviour depends 
crucially on whether the wavelength (A) is greater than 
the average distance between air blocks. When (b)/X is 
less than one, the localization effect is weak. We also 
observe that with the added disorder, the transmission 
is enhanced in the middle of the gaps. Similar enhance- 
ment due to disorder has also been reported recently jl9) . 
Differing from fl^| , however, the transmission at frequen- 
cies within the gaps of the corresponding periodic arrays 
in the present system is not always enhanced by disor- 
der. Instead the transmission is reduced further by the 
disorder near the band edges. 
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FIG. 3. Diagrams (a), (b), and (c) show the Lyapunov ex- 
ponent (LE) in broken lines and its variance in solid lines as 
a function of k{b)/n. Diagrams (d), (e), and (f) present the 
plots of the exponent versus its variance. Here (3 = 10 -4 . 

Fig. [| presents the results for Lyapunov exponent (LE) 
and its variance as a function of non-dimensional fre- 
quency k{b)/ir. At low disorders, the variance of LE in- 
side the gaps is small. Contrast to the optical case [ p"7| , 
there are no double maxima inside the gap. With in- 
creasing disorder, double peaks appears inside the allow 
bands. When exceeding a certain critical value, however, 
the double peaks emerge. The higher frequency, the lower 
is the critical value. For example, the double peaks are 
still visible in the first allow band (c. f. Fig. ||(b)), while 
there is only one peak inside the higher passing bands. 



Meanwhile, the increasing disorder reduces the band gap 
effect and smears LE, in accordance with Fig. |[ We 
also plot LE- variance relation in Fig. ||. With increasing 
disorder, we do not observe genuine linear dependence 
between LE and its variance, as expected from the single 
parameter scaling theory. 

In the past the localization phenomenon in ID is usu- 
ally characterized by LE, here we propose to use a phase 
behavior of the waves to characterize the localization. We 
compute the energy density inside the sample from 
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where m refers to either water or air. Meanwhile, the 
phases of p(x) and u(x) are recorded. Expressing p(x) 
and u(x) as A p (x)e l6p ^ and A u {x)e l6ulyX \ we construct 
unit phase vectors v p = cos6 p e x + sm.9 p e y and v u — 
cos9 u e x + sinf^ey. Then the behavior of the phase vec- 
tors along the path are investigated. Physically, these 
phase vectors represent the oscillation behavior of the 
system. 

Typical results for the spatial distribution of the en- 
ergy and the phase behaviour for the given disorder and 
air fraction are shown in Fig. ^. Here for the sake of con- 
venience, only the phase vectors at the interfaces between 
air and water are shown. First, we note that the energy 
density is constant in each individual block. This is a 
special feature of ID classical systems, and can be veri- 
fied by a deduction from Eq. (|^) . We find that when the 
sample size is sufficiently large, waves are always localized 
for any given amount of randomness. When localized, the 
waves are trapped inside the medium, but not necessarily 
confined at the site of the source, unless the band gap ef- 
fect is dominant. The energy distribution does not follow 
an exponential decay along the path. This differs from 
situations in higher dimensions [||,^0). It is also shown 
that the energy stored in the medium can be tremendous 
(c. f. Figs. [|(a) and (b)). With increasing sample size, 
the peak amplitude may grow, pointing to the stochastic 
resonance, in agreement with When disorder is weak 
and for frequencies within the gaps, such a stochastic 
resonance behavior disappears. 

It is also observed that for all frequencies, there is a 
collective behavior for the phase vectors. In Fig. |J(c) we 
plot the phase vectors in three spatial regimes, namely, 
near the transmitting site, in the middle of the duct, and 
at the far end from the source. Symbols p L , p R , u L , u R 
appearing in Fig. ||(c) denote respectively the phase vec- 
tors for the pressure and the velocity fields on the left and 
right side of the air blocks. It is clear that when waves are 
localized, all the phase vectors of the pressure field are 
pointing to either 7r/2 or —ir/2, and perpendicular to the 
phase vectors of the velocity field. The pressure at the 
two sides of any air block varies in phase. Mostly, the 
two sides also oscillate in phase. Different from higher 
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dimensional cases in which all phase vectors of localized 
fields point to the same direction ||[20) , the present phase 
vectors are constant by domains; this ensures no energy 
flows. The velocity field in neighboring domains oscil- 
late exactly out of phase. The phase vector domains are 
sensitive to the arrangement of the air blocks. We stress 
that such a phase ordering not only exists for the bound- 
aries of the air blocks, but also appears inside the whole 
medium. 



of such a collective phenomenon may be regarded as an 
indication of a kind of Goldstone modes in the context 
of the field theory Q. 
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FIG. 4. Energy distribution and phase diagrams for local- 
ized state. Here (3 = 10 -4 , A = 0.3 and kip) = 0.3-7T is a 
frequency in the first allowed band, (a) Energy densities in 
the air blocks, (b) Energy densities in water, (c) Phase vec- 
tors at the interfaces for three spatial ranges of the medium. 
The unit of energy density is J/m 3 . 

The coherence behavior is a unique feature for wave 
localization, which could be verified by the cross correla- 
tion measurement. At the far end of the sample, however, 
the phase vectors become gradually disoriented, implying 
that the energy can leak out only at the boundary due 
to the finite sample size. This boundary effect vanishes 
exponentially as the sample size is increased. The fact 
that the phase vectors are constant in domains indicates 
that once waves are localized, no more energy can be 
pumped into the system. When localization is evident, 
increasing the sample size by adding more air blocks will 
not change the patterns of the energy distribution and 
phase vectors. Therefore the energy localization and the 
phase behavior are not caused by the boundary effect. 

In summary, we have demonstrated a phase transition 
in acoustic propagation in an ID random liquid medium. 
The results have shown that waves are always confined 
in a finite spatial region. The disorder leads a significant 
energy storage in the system. It is also indicated that 
the wave localization is related to a collective behavior 
of the system in the presence of multiple scattering, also 
observed for higher dimensions pj2Q]. The appearance 
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